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Abstract. We consider quantum spin systems defined on finite sets V equipped witli a metric. In 
typical examples, V is a. large, but finite subset of Z**. For finite range Hamiltonians with uniformly 
bounded interaction terms and a unique, gapped ground state, we demonstrate a locality property 
of the corresponding ground state projector. In such systems, this ground state projector can be 
approximated by the product of observables with quantifiable supports. In fact, given any subset 
X C V the ground state projector can be approximated by the product of two projections, one 
supported on X and one supported on X'^, and a bounded observable supported on a boundary 
region in such a way that as the boundary region increases, the approximation becomes better. Such 
an approximation was useful in proving an area law in one dimension, and this result corresponds 
to a multi-dimensional analogue. 

1. Introduction 

The intense interest in entangled states for purposes of quantum information and computation 
during the past decade has stimulated new investigations in the structure of ground states of 
quantum spin systems. It was soon found that the focus on entanglement, and the mathematical 
structures related to this concept, also provide a new and useful way to investigate important 
physical properties relevant for condensed matter physics [HHO]. In particular, progress was made 
in our understanding of constructive approximations of the ground states of quantum spin models, 
and about the computational efficiency of numerical algorithms to compute ground state properties 
and to simulate the time evolution of such systems |18^ [T9] . 

It is now understood that there is a relationship between the amount of entanglement in a state 
and the degree of difficulty or the amount of resources needed to construct a good approximation 
of it [T7]- To make this more precise one has to use a quantitative measure of entanglement. For 
pure states the entanglement entropy is a very natural such measure for bi-partite entanglement. 
For a two-component system with components A and B (e.g., A and B label a partition of the spin 
variables into two sets), the entanglement entropy of the system in a pure state tp with respect to 
this partition is defined as the von Neumann entropy of the density matrix describing the restriction 
of "0 to subsystem A. Mathematically, this density matrix, pA, is obtained as the partial trace of 
the orthogonal projection onto the state ip calculated by tracing out the degrees of freedom in B. 

It was conjectured that under rather general assumptions, the entanglement entropy associated 
with A satisfies an Area Law. By Area Law one means that if A is associated with a subvolume 
of a fixed system in a ground state then, with a suitable definition of boundary, the entanglement 
entropy of A should be bounded by a constant times the size of the boundary of A (e.g., the surface 
area in the case of a three-dimensional volume). 

One way to understand this conjecture is to observe that it holds almost trivially for the so-called 
Matrix Product States (MPS), aka Finitely Correlated States [3] and their generalizations in higher 
dimensions known as PEPS [16]. It is therefore plausible that to prove the Area Law it is key to 
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have a good handle on MPS-hke approximations of a general ground state. By using such a strategy 
Hastings obtained a proof of the Area Law for one-dimensional systems [5] . 

The goal of this paper is to prove a generalization to arbitrary dimensions of the approximation 
result Hastings used in his proof of the Area Law in one dimension. This is Theorem 12.21 below, 
which demonstrates that the ground state projector of certain gapped quantum spin systems can 
be approximated by a product of three projections with known supports. We describe this result 
in detail in the next section. 



2. A Ground State Approximation Theorem 

We begin with a brief description of the quantum spin systems that will be considered in this 
work. Let V be a countable, locally finite set equipped with a metric d. In most examples, V = 
and d is, for example, the Euclidean metric. Since our arguments do not make any use of the 
structure of the underlying lattice, we present our models in this more general setting. Most of 
our analysis applies on finite subsets V C V. More specifically, to each x £ V, we will associate a 
finite-dimensional Hilbert space TCx', the dimension of TCx will be denoted by n^. Set to be 
the complex Ux x Ux matrices defined over TCx- We denote by TCy = 'S)xgv'^^ Hilbert space 
of states over V, and similarly Ay = 'S)xgV ^^i^ algebra of local observables. For any finite 

subsets X d Y <Z V , the local algebra Ax can be embedded in Ay by using that any A S Ax 
corresponds to yl I S Ay- Thus, the algebra of quasi-local observables, .Ay, can be defined as 
the norm-closure of the union of all local algebras; the union taken over all finite subsets of V. 
Moreover, we say that an observable A G Av is supported in X C V, denoted by supp(^) = X, if 
X is the minimal set for which A can be written ss A = A®1 with A G Ax- 

A quantum spin model is then defined by its interaction and the corresponding local Hamilto- 
nians. An interaction is function $ from the set of finite subsets of V into Ay with the property 
that given any finite y C V, ^iy) = ^(V)* G Ay- Given such an interaction, one can associate a 
family of local Hamiltonians, parametrized by the finite subsets of V, defined by setting 

(2.1) Hv=Y. HX) 

xcv 

for finite V C V- Since each local Hamiltonian Hy is self-adjoint, it generates a one parameter group 
of automorphisms, which we will denote by tY , that is often called the finite volume dynamics. 
This dynamics is defined by setting 

(2.2) (A) = e''^'' Ae-''^'' , AeAy- 

When the subset V on which the dynamics is defined has been fixed, we will often write Tt to ease 
the notation. 

Theorem 12 . 2 1 b elow holds for certain quantum spin models. Our first assumption pertains to the 
underlying set V. We must assume that there exists a number /io > for which 

(2.3) K^, = sup V e-^«"'(^'^) < oo . 
Next, we make the assumptions on the interactions we consider precise. 

Al: We consider interactions that are of finite range. Specifically, we assume that there exists 
a number R > 0, called the range of the interactions, for which ^{X) = if the diameter of X 
exceeds R- Here, for any finite X C V, the diameter of X is 

(2.4) diam(X) = max{d{x, y) : x,y £ X} 
where d is the metric on V. 

A2: We will only consider uniformly bounded interactions. Hence, there is a number J > such 
that for any finite X C V, \\<^>{X)\\ < J- 
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A3: Let be the characteristic function defined over the set of finite subsets of V, i.e. for any 
finite X C V, x^i^) = 1 if ^{X) ^ 0, and x^i^) = otherwise. We assume the following quantity 
is finite 

(2.5) iV$ = sup ^ |X|x<i.(X) < cx) 

where the sum above is taken over finite subsets X G V and \X\ is the cardinality of X. Another 
relevant quantity, which often appears in our estimates, is 

(2.6) C$ = sup y x<fiX) < A$ . 

xex 

A4: We will assume that a given local Hamiltonian Hy has a unique, normalized ground state 
which we denote by tpQ. Pq will be the projection onto this ground state in TCy, and we will label 
by 7 the length of the gap to the first excited state. 

Much progress has been made recently in proving locality estimates for general quantum spin 
systems, see e.g. [6l[lT]. The following estimate was proven in |14j . 

Theorem 2.1 (Lieb-Robinson Bound). Let V be a countable, locally finite set equipped with a 
metric d for which 1^2. 3\) holds. Let $ be an interaction on V that satisfies assumptions Al, A2, 
and A3. For every /i > /io, there exists numbers c and v such that given any finite V C V and any 
observables A G Ax and B G Ay with X,Y C V and X CiY = f}>, the bound 

(2.7) ||[ri^(A),B]|| <c Pillion min[|9$X|,|9$y|]e-'^(^(^'^)-^'l*l) 
holds for all t gM. 

Here for any finite X C V, the set 

(2.8) d^X = {xeX : there exists Y CV with x €Y,Y nV\X ^0, and, <^{Y) / 0} , 

the number d{X, Y) = mm{d{x, y) : x £ X,y £ Y}, and we stress that both c and v are independent 
of F C V. 

Before we state Theorem 12. 2( some further notation is necessary. Let ^ be an interaction which 
satisfies assumptions Al - A4 above. Again, most of our analysis applies on finite sets V C V, and 
often the particular finite set under consideration is the one whose existence is guaranteed by A4. 
In general, for finite sets X C V C V, we will denote by dX the R-boundary of X in V, i.e., 

(2.9) dX = {x £ X : there exists y £V\X, with d{x, y) <R}. 

We note that with X fixed, the set dX is independent of V, for V sufficiently large, and so we 
suppress this in our notation. To state our main result, we fix an additional subset X d V . The 
following sets, defined relative to this fixed set X and dependent on a length-scale £ > R, will also 
play an important role. Set 



(2.10) Xint = X^l^{£) = {x£X :iov all y G dX, d{x, y) > £} , 

(2.11) Xbd = X^^{i) = {x£V : there exists y G dX, d{x, y) < i} , 
and 

(2.12) X^^t = '^cxt W = {x G \ ;f : for all y G dX, d{x, y) > £} . 



The sets Xi^t, X\,^, and Xc^t correspond respectively to the -^-interior, the ^-border, and the 
^-exterior of in y. We note that for any finite V C V and each X C V, the sets Xi^t, Xi,^, and 
^ext are disjoint, and moreover, V = Xi^t U X^^ U X^xt- We caution that the notation dX^xt here 
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refers to the boundary of this set as a subset of V, as indicated above; not as a subset of V. It will 
also be important for us that there exists a number C > 0, independent of both i and V, for which 

(2.13) max{|a^int|,|^bd|,|5;fext|} < Ci\dX\. 

The above inequality constitutes the main structural assumption on the set V and the reference set 
X C V. We may now state our main result. 

Theorem 2.2. Let V be a countable, locally finite set equipped with a metric d for which 12. 3\) holds, 
and moreover, let $ be an interaction on V which satisfies assumptions Al - A4, take V C V to 
be the finite set from A4, and suppose X d V satisfies i2.13\) . For any i > R, there exists two 
projections Px G Ax and Px<= G A.v\x and an observable Px^d ^ •^A'bd{3^) ^^^^ ll-^'^bdll — ^ such 
that 

(2.14) WPx^^PxPx^ - PoW < KC^N^\dX\^f/^e-'/^^. 
Here ^ is defined by 

for /X > 2/io. 

In many situations, it is useful to have a non-negative approximation of the ground state pro- 
jector. For this reason, we also state the following corollary. 

Corollary 2.3. Under the assumptions of Theorem \2. 21 the estimate 

(2.16) WPx^PxPXdPx^dPxPx^ - P4 < KC^m\dX\^l"^e~''^^ 

also holds. 

Since all the observables involved have norm bounded by 1, Corollarv 12.31 is an immediate con- 
sequence of Theorem 12.21 

With the methods discussed in [H], Theorem 12.21 can easily be generalized to other types of 
interactions. For example, the finite range condition Al can be replaced with an exponentially 
decaying analogue. Moreover, the uniform bound on the local interactions terms, i.e. A2, can also 
be lifted. One could also consider systems where the underlying sets satisfied a bound of the form 
(j2.13p with a larger power of L This would increase the power of £ in the statement of Theorem l2.21 
but not effect the exponential decay. We give the argument in the context described above for 
convenience of presentation. It would be very interesting to see that the methods developed here 
can be used to prove an area law, similar to Hastings' result [5], in dimensions greater than 1, but 
this is beyond the scope of the present work. 

3. Basic Set-Up 

Before beginning the proof of this theorem, we will introduce the basic objects to be analyzed. 
Let y C V be the finite set described in A4, take £ > R, and fix X C V satisfying (j2.13p . Note 
that Hy, as defined in (12. Ih . can be written as a sum of three local Hamiltonians 



(3-1) Hv = H\.^^ + Hx^d + ^A'cxt ' 

where Hx^^ ™ (j2.ip with V = A'bd, and for Z G {^^int, '^cxt}) 
(3.2) E 

XCV: 

corresponds to local Hamiltonians that include boundary terms. Without loss of generality, we may 
subtract a constant, namely {^q, Hyipo), and thereby assume that the ground state energy of each 
of these Hamiltonians is zero. In this case, each term in (|3.ip has zero ground state expectation. It 
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is not clear, however, that each of these terms, when apphed to the ground state, have a quantifiably 
smah norm. To achieve smaU norms, we introduce non-local versions of these observables. 

Generally, given any self-adjoint Hamiltonian H on Tiy, a local observable A E Ay, and a > 0, 
we may define a new observable 

(3.3) (A)^ = V?/" 

where here Tt{A) = e**^yle~**^. It is easy to see that for every a > 0, ||(A)q,|| < ||^||. 

In the present context, we will define the (•)„ operation with respect to the Hamiltonian Hy 

rY 

Hamiltonian, i.e., 



and its corresponding dynamics tY ■ Observe that these non-local observables still sum to the total 



(3.4) Hy = {H%J^ + {Hx,X + (^Lt)- 
They also have zero ground state expectation, i.e., 

(3.5) (^o,(^L)"^o) = (^o,(^;t,Ja^o) = (^o,(^LJ"^o) = 0. 

We will prove, in Proposition 15.21 that each of these smeared-out Hamiltonians, when applied to 
the ground state, has a norm bounded by an exponentially decaying quantity. The cost of this small 
norm estimate is a loss of locality, which is expressed in terms of the support of these Hamiltonians. 

Equipped with the Lieb-Robinson bound in Theorem 12. H we begin the next step in the proof of 
Theorem 12.21 Here we approximate each of (^T^.^Jcn {Hx^i)a, and (iJ^^^J^ with local observables 
Mx{ol), Mx^^(2i){o') and Mx<:{a) respectively. These observables will be constructed in such a way 
that they not only have a quantifiable support but also a small vector norm, when applied to the 
ground state. 

To make these approximating Hamiltonians explicit, we introduce local evolutions. Define three 
different dynamics, each acting on Ay, by setting 

Tt^{A) = e'^^^Ae-'*^^ , 

(3.6) Tf'"'^^'^\A) = e**-^-*bd(2'^)^e"**-^*bd(2^) , 

Tf"{A) = e'^Hv\x Ae-^*^y\x , 

for any local observable A S Ay and any t G M. The set X\^^[2€) is as defined in (j2.1ip with the 
length scale doubled, and the local Hamiltonians used in the evolutions above are as in (j2.ip . The 
approximating Hamiltonians are given by 

MMa) = ^JT^rr {H^J e"^^' dt . 

The observables defined above, i.e. in (|3.7|) . have been chosen such that supp(M;i'(q)) = X, 
supp(M;^;'j^^(2£)(ct)) = '^bd(2^), and supp(M;t:'c(a)) = V \ X . To ease notation, we will often write 
^x^ad^Oi) = M p^^^^2t){o') and suppress the exact size of the support. 

The following technical estimate summarizes the results mentioned above. 

Lemma 3.1. Let V he a countable, locally finite set equipped with a metric d for which 112. holds, 
and let $ be an interaction on V which satisfies assumptions Al - A4. For the finite set V G V 
from A4 and any set X C V satisfying 112.1 3\) for all £ > R, the estimate 

(3.8) \\Hy - {Mx{a) + Mx^^ia) + Mxcia))\\ < K\dX\i^/^e~l , 
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holds along the parametrization 2ae£ = fiv'^ where fi > 2fiQ. The numbers ^ and e are defined in 
terms of the gap 7 and the quantities /i and v from the Lieb-Robinson estimates as 

(3.9) < - = (1 - e)^ = ^/i . 

Along the same parametrization, the bound 

(3.10) max{||M;f (a)^oll , \\M;,^A»)M , ||M;fc(a)V'o||} < K\dX\f/^e-l. 
also holds. 

These estimates play an important role in the proof of Theorem 12.21 We prove this lemma in 
Section [5l For now, we use it to prove Theorem! 



4. Proof of Theorem 12.21 



The first step in verifying the bound claimed in (j2.14p is to find an explicit approximation to the 
ground state projector Pq. With this in mind, define 



(4.1) Va = 

for any a > 0. Clearly, Va G -^v and for any vectors f,g € Tiy, the spectral theorem implies 



(/, (Va - Po)g) 



e'^'d{f,Exg)dt - {f,Pog) 



dtd{f,Exg) 



(4.2) 



e 40, d{f,Exg), 



where we have denoted by Ex the spectral projection corresponding to Hy- This readily yields 
that 



(4.3) 



\Va-Po\\ < e-^, 



where 7 is the gap of the Hamiltonian Hy- 

Using the operators introduced in (j3.7p . we define an analogous ground state approximate by 
setting 



(4.4) 



Here we have dropped the dependence of Mx^, ^x^^, and Mx'^ on a. Clearly, 



(4.5) 



Va-Po 



< 



+ \\Va-Po\\. 



The final term above we have bounded in ()4.3p . To bound the first term in ()4.5p . we introduce the 
function 



(4.6) 

One easily calculates that 
(4.7) F/(A) = -it 



Ft{X) = e^HMx+Mx^^^+Mxc)t ^i(l-X)Hvt_ 



jxiMx+Mx^^+Mxc)t s^Hy - (Mx + Mx^^ + Mxc)} e^(i-^)^^*. 
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Using Lemma l3.ll we conclude that, if latl = fiv'^, then 

\\Ft{l) - Ftm\ e-^'"dt 

^ J \\Hv - {Mx + Mx^^ + Mx^ 
(4.8) < K\dX\fe~'i . 







T^a 'Pa 





< 



\t\ e 



dt 



In the bound above, and for the rest of this section, the number K which appears in our estimates 
will depend on the parameters of the quantum spin models, but not on i. It may change from line 
to line, but we do not indicate this in our notation. 
Inserting both (fOj) and ([i^ into yields 



(4.9) 



again along the parametrization 2aei = fj.v'^. 
Set 



(4.10) 



6 = ^3/2 



e 2«, 



The projections Px and Px<:, which appear in the statement of Theorem 12.21 defined to be the 
spectral projection corresponding to the matrix Mx, respectively Mx<:, onto those eigenvalues less 
than 5. By definition, Px G Ax and Px'^ G Ax'^ as claimed. Moreover, using Lemma |3. 11 it is easy 
to see that 



(4.11) 



max \\(1 - Pz) ipoW < max - ||Mz'0oll < KldXle 2« . 
ze{x,X''} z&{x,x<'}5 



by construction; hence this choice of 6. 

With this bound, we can insert these projections into our previous estimates. In fact. 



(4.12) 



PaPxPx'^ — Pq 



< 



Va - Po]PxPx 



+ ||Po (1 - PxPx^ 



The first term above is estimated using (14. 9p above. Since 

1 



(1 - PxPx^) = 2 - ^'^)(^ + Px^) + (1- Px^)il + Px)} , 



(4.13) 
it is clear that 

(4.14) ||Po(l - PxPx^)\\ < \\Po{l-Px)\\ + \\Po{l - Px^)\\ < 2K\dX\e-'/^^, 

using ()4.1ip . Therefore, we now have that 



(4.15) 



VaPxPx^ 



again along the parametrization 2aei = ^v^ . 

In order to define the final observable Pxy,^, we write 



(4.16) 



VaPxPx^ 



^i{Mx+Mx +Mxc)t -i(Mx+Mxc)t i{Mx+Mx<^)t 



PxPx^e 



dt. 



Since the supports of Mx and Mx<^ are disjoint, it is clear that 



J{Mx+Mxc)t 



(4.17) 



PxPx^ - PxPx^ 



l{e'^^^Px - PxW'^'^'Px^ + Px^) + \{e'^'''^'Px^ - Px^W'^'Px + Px) • 
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Moreover, for Z G we have that 



\{f,ie 

(4.18) 

and therefore, 
(4.19) 

If we now define the operator 
(4.20) Va 



iMzt 



^)Pzg)\ 





< 5\t 



< 26\t\ . 



J{Mx+Mx^^+Mx<^)t^-i{Mx+Mxc)t^-at^ 



then we have just demonstrated that 

\\Pc.PxPx^ -P4 < WPaPxPx^ - APxPx4 + WVaPxPx^ " ^ol 



< 26 J ^ 

TT 



t|e-"* dt + K\dX\£^e~ 



2€ 



(4.21) 



< K\dX\fe~^. 



The proof of Theorem 1 2. 2 1 is complete if we show that the operator Va, defined in (j4.20p above, is 
well approximated by a local observable. In fact, below we introduce an observable -P^'bd ^ •^Xi,^{3e) 
with ll-P^'bdll — ^ which 



(4.22) 



^bd 



Pn 



< K\dX\''el^e 



From this estimate. Theorem 12.21 easilv follows as 



\Px^^PxPx^ - P4 < 



Px^, - Va PxPx 



+ WVaPxPx^ - Po 



(4.23) 



We need only prove that such an observable Pxi,d exists. 

The existence of this observable is simple. We just take the normalized partial trace of Pa over 
the complimentary Hilbert space, i.e., the one associated with ^^^(S^)'^ C V. To prove the estimate 
in ()4.22p . it is convenient to calculate this partial trace as an integral over the group of unitaries, 
technical details on this may be found in [2] and also [13]. We recall that given an observable 
A G Ay and a set Y C V, one may define an observable 



(4.24) 



{A)y 



U*AUfi{dU), 



where U{Y^) denotes the group of unitary operators over the Hilbert space Wye and fi is the 
associated, normalized Haar measure. It is easy to see that for any A £ Ay, the observable {A)y 
has been localized to Y in the sense of supports, i.e., {A)y G Ay- 
Let us define 



(4.25) 



{Pa)Xi,ii3e)- 

Clearly Px^^^ G ^A;bd(3£) and \\Px^,J\ < 1 as desired. 

The difference between Va and Pf^d ^^n be expressed in terms of a commutator, in fact. 



(4.26) 



To 



Px. 







[ u* 


Va,U 


Ju(Xi,i{3er) 





fi{dU). 
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Thus, to show that the difference between Va and Px^^ is smah (in norm), we need only estimate 
the commutator of Va with an arbitrary unitary supported in Abd(3^)'^. 
This follows from a Lieb-Robinson estimate. Note that 



(4.27) 



dt. 



To estimate the integrand, we define the function 



(4.28) 



fit) 



J{Mx +M;t-bd +Mxc)t^-i{Mx+Mxc)t jj 



A short calculation demonstrates that 



(4.29) 



nt) 



^i(Mx+Mx^^+Mxc)t j^^^^^-i{Mx+Mxc)t ^ jj 

The form of the derivative appearing in (I4.29P suggests that we define the evolution 

(4.30) at{A) = em^+Mx^^+Mxc)tjs^^-mx+Mx^^+Mxc)t ^ ^ 

With this in mind, we rewrite 

f(t) = i (^{'^'ix+Mx^^+Mx<^)tyj^^^^^-i(,Mx+Mxc)t 

= i Qt(M;,,Je^(*^^+^^bd+*^^^)*e-^(*'^^+*^^=)*,?7" 

(4.31) = iat{Mx^,)f{t) + i [at{Mx^,), U] e^(^^^+A^^bd+*^^=)*e-*(*^^+^^^)*. 

Written as above, the function / can be bounded using norm-preservation, see [TT] and [H] for 
details. For example, let V{t) be the unitary evolution that satisfies the time-dependent differential 
equation 

d 



ijV{t) = V{t)at{Mx^^) with V(0) = 1. 



^i{Mx+Mxc)t ^-i{Mx+Mx^^+Mx<^)t _ 



(4.32) 

Explicitly, one has that 
(4.33) V{t) 
Considering now the product 
(4.34) 
it is easy to see that 

(4.35) g'{s) = V'{s)f{s) + V{s)f'{s) = iV{s) [a,(M^,J, U] V{sY. 

Thus 



9{s) = V{s)f{s) 



(4.36) 

and therefore, 
(4.37) 

The bound 
(4.38) 

readily follows. 



V{t)f{t)=g{t)-g{^) 



g'{s)ds, 



fit) = iVity / V{s)[asiMx^,),U]Visrds. 



\\f{t)\\< f'\\[as{Mx,,),U]\\ ds. 
Jo 
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Unfortunately, the Lieb- Robinson velocity associated to the dynamics at{-) grows with £. For this 
reason, we estimate ()4.38p by comparing back to the original dynamics. Consider the interpolating 
dynamics 



(4.39) 



hs{r) = Ur {Ts-r{A)) 



for any local observable A and < r < s. Here Tj-^A) = e^^^^Ae ^^vr\ With s fixed, it is easy to 
calculate 

(4.40) h',{r) = ias-^r {[{Mx + M^^^ + M;,.) - Hy, MA)]) . 
We conclude then that 

\\as{Mx^,)-Ts{Mx^,)\\ = [ K{r)dr 

Jo 

< 2 \\Mx^J \\Hv - {Mx + Mx^, + Mx^)\\ s 

< 2(JC$|A'bd|) {K\dX\l^'^e~~i^ s 

(4.41) < K\dX\'^l^'^e~~i s, 
where, for the second inequality above, we used Lemma 13.11 and the bound 



< 



< J C(j)| A^bdl ■ 



(4.42) 

From (j4.38p . it is clear that 

(4.43) ||/(t)||< r \\\rs{Mx,,),U]\\ ds + 



rf- ($(X)) 



dt 



[as{M, 



The first term above, we bound, using Theorem 12. H as follows 



{T,{Mx^:)m\ds < C\dXU2i)\\\Mx^ 



1*1 



Ts{Mx,,),U]\\ ds. 



-tiidiX^a{2e),U)-v\s\) 



(4.44) 



1 



For the second, we apply ()4.4ip . Thus, we have found that 



\\[Va,U]\\ < \ - 



e-'^' dt 



t^e-"*' dt 



fJ,V \ TT 



(4.45) < K\\dX\^f/'^e~~i . 

Combining this with (j4.26p . we obtain that 



(4.46) 



and the proof of Theorem 12.21 is now complete. 
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5. A Technical Estimate 
In this section, we will prove Lemma |3. 11 For convenience, we restate the result. 

Lemma 13. 11 Let V be a countable, locally finite set equipped with a metric d for which \2. 3|) holds, 
and let $ be an interaction on V which satisfies assumptions Al - A4. For the finite set V <ZV 
from A4 and any set X <ZV satisfying i2.13\) for all i > R, the estimate 



(5.1) 



\Hv - (Mxia) + Mx^^ia) + Mxc{a))\\ < K\dX\£^/^e- 



holds along the parametrization 2ae£ = fiv'^ when > The numbers ^ and e are defined in 

terms of the gap 7 and the quantities fi and v from the Lieb-Robinson estimates as 



(5.2) 



0<- = (l-e)/i 



7 

09, 9 A* • 

^ZyZ _|_ ryZ 



Our bounds on the prefactor K depend on various parameters, e.g. J'^ ,C^, N^, fj,,v,j, R, but 
it is independent of i,V, and X. Along the same parametrization, the bound 



(5.3) 

also holds. 



max{||M;,(a)^o||,||M;,,,(a)Vo||,||M;,c(a)V'o||} < K\dX\el^e—^. 



To prove Lemma 13.11 we begin with a few simple propositions. Let V <Z V he the finite set 
described in A4, take £ > R, and fix A' C ^ satisfying (j2.13p . As we demonstrates in Section O 
see also Section [21 the local Hamiltonian corresponding to <I> in y can be written as 

(5-4) Hy = H';,^^^ + Hx^^ + H'x^^^ 

where the sets ^mt, '^bdj and X^xt, and the corresponding local Hamiltonians, each depend on a 
length scale i > R. We begin with a basic commutator estimate. 

Proposition 5.1. Let $ be an interaction on V which satisfies assumptions Al - A4. Let V C V 

be the finite set described in A4, take i > R, and fix X C V satisfying i2.13\) . One has that 

(5.5) max[\\[Hv,H'^J\\A\[Hv,H;,J\\A\[Hv,H';,^J\\}< 

where the quantities C$ and are as introduced in (j2.6p and (j2.5p . 
Proof. Define 



(5.6) 



int 



Cbd — [HviHx^^] , and C, 



cxt 



From the definitions, (|2.10p . (|2.1ip . and (|2.12p . and the fact that ^ > i?, it is clear that Qnt 
[Hx^^,H^x,J^ C-cxt = [Hx^^,H^x,J, and -Cbd = Qnt + Cext- Clearly, 



(5.7) 

and therefore, 



Cn 



XCV: 



int I 



< 



E 



< 



(5.8) 



XGdXint XCV: 

E EEEii[^(n^>w]i 

x€dXint XCV: Z£X YCV: 

^ E EEE^*(^)^*w 

< 2j2C$iV$|5A^int| 



XG^A^int XCV: Z&X YCV: 

xGX zeY 
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A similar estimates applies to Ccxt- Using (j2.13p . the proof is complete. □ 

As is discussed in Section O the first step in the proof of Theorem [22] is to introduce the smearing 
operation (•)q,, see ()3.3p . One consequence of this definition is 

(5.9) Hv = {H\J^ + {Hx,X + (^LJ- 

When applied to the ground state, each of the terms above can be estimated, in norm, in terms of 
the system's gap. This is the content of Proposition 15.21 below. 

Proposition 5.2. Let ^ he an interaction on V which satisfies assumptions Al - A4. For any 

finite set X CV (ZV, a> and A £ {H^,^^^, H^i,^, H^^^J , 

(5.10) ||(^)„^o|| < &^e-2^. 

7 

Proof. Let A E {H^, , H^^^, H^^^^}. As previously discussed, see (jS.Sp . = {ipo, AiI^q) = {tpQ, (j4)a^o)- 
Setting Ca_ = [Hy ,A], we find that 

\\{A)aM < -\\Hv{A)aH\ 

7 

= -\\[Hv,iAUM 

7 

(5.11) = -\\iCA)aM- 

7 

The value of ||(CA)aV'o|| can be estimated using the spectral theorem. For any vector /, one has 
that 

(/,(CA)aVo) = r e"''\f,T^iCA)i^o)dt 



oo 



-oo 
oo 



TT 

a 
vr 

oo r~~ /-oo 



e-"* / e''^d{f,ExCA^o)dt 

oo ./ 7 

^ P e-'^'^e^'^dtdU^E^CAi^o) 



(5.12) = / e-*^d{f,E^CA^(^). 

In the third equality above we have introduced the notation E\ for the spectral projection corre- 
sponding to the self-adjoint operator Hy, and we also used that {ipQ, CAi^o) = 0. The last equality 
is a basic result concerning Fourier transforms of gaussians (and re-scaling), (see e.g. (5. 59) in |13j). 
We conclude then that 

(5.13) \{f,iCA)aM < e-^||/||||C7^||, 
and with / = (Ca)«V'o we find that 

(5.14) KCaUMI < e~^\\CA\\. 
Putting everything together, we have shown that 

-- 

(5.15) ||(^)a^0|| < -\\iCA)aM < 

7 7 

□ 
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The next step in the proof of Theorem 12.21 is to make a strictly local approximation of the 
smeared terms appearing in (j5.9p above. We defined these local approximations in Section [3l see 
(j3.7p . and Proposition 15.31 below provides an explicit estimate. 



Proposition 5.3. Let V he a countable, locally finite set equipped with a metric d for which i2. 3\) 
holds, and let be an interaction on V which satisfies assumptions Al - A4. Let V C V be the 
finite set described in A4, take i > R, and fix X C V satisfying \2.13\) . The estimate 



max ■ 



(5.16) 



< 2C- 



-\dX\ 4 



2e/i 
vr 



^3/2 



3^lR 



e i. 



holds along the parametrization 2aei = ^v^ where fi > 2/io. Here 



(5.17) 



< - = (1 - e)/i = — 2 M ■ 



Proof. We will prove the estimate for (i^^. — Mx{oi), the other bounds follow similarly. Note 
that 



(5.18) 



{H\J. - Mx{a)\ < ^ I^JrtiH'xJ " ^,^(^1 



-at' 



dt. 



To bound the integral above, we introduce a parameter T > 0. For |t| > T, we use that 



n H. 



< 



t\ 



ds 



(5.19) 

From this, it follows readily that 

t\>T 



< n[Hv,H%J\\\t\ 



(5.20) 



where we used Proposition 15.11 

For |t| < T, the estimate below is an immediate consequence of Lemma 3.1 in 

\t\ 



(5.21) 



< 



ds. 



The above commutator may be written as 



(5.22) 



E 



Y. [^(xj.Tf (4.(y'))] 



XCV: Y<ZV: 



Estimating this, we find that 



Hy — Hx,T. 
(5.23) 



H 



Xi„ 



x€dX J/eA'int XCV: YCV: 

xex yeY 



^ ^"^^ E E E E 



■^l{d{X,Y)-v\s\) 



xedX yeXint XCV: YCV: 

xGX y& 



where we have used the Lieb- Robinson bound, i.e. Theorem 12.11 with /i > /io from (j2.3p . Since ^> 
has a finite range i?, it is clear that d{x,y) — 2R < d{X,Y). This implies that 

^x,^,j]\\<cJ'CMe"''' E e-'^('^(^'^)-''l^l) . 

xf^dX yeXint 



(5.24) 



Hy - LLxiT^ ( H. 
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Now, for each fixed x G dX, d{x, y) > i — 2R. Summing on all y yields 



(5.25) 



Hy — Hx,T^ ( H 



< 



cK^/2J'^C<s>N<s>e^>'^\dX\e 



-f,{i/2-v\s\) 



for fi > 2fiQ. Comparing back to (j5.2ip . this bound demonstrates that 



(5.26) 



^tj.v\t\ _ -^^ 
jJLV 



and hence, for any T > 0, we have that 



a 



(5.27) 



< 2cK^/2e' 



/if 



\dX\e 



Adding our results for large and small T, it is clear that 



(5.28) 



< 



+ CK^i^e-^^ ■ exp 



~2a 



We now choose a parametrization which illustrates the decay. Let a satisfy 



~2a 



el with < e 
and choose T to satisfy the equation 
(5.29) 

In this case, 2aT^ = (1 — e)/z£, and thus 
(5.30) 



1 + 



7 



< 1, 



[IV 

~2oL 



0. 



fiv \ V vr 



^3/2 + 



,3fiR 



Equipped with theses estimates, we are now ready to prove Lemma l3.ll 

Proof, (of Lemma 13. II ) Using (j5.9p and Proposition 15.31 we find that 
II - {Mx{a) + Mx^^ia) + M;t=(a))|| 

(5.31) 



< 6C 



/if \ V vr 



7 



-e 4q 



(5.33) 



7 



e « 



along the parametrization 2a€i = /if^. This is the bound claimed in (j5.ip 
To see that ()5.3p is true, note that, for example, 

(5.32) ||M^(a)V'o|| < UH^xJaM + II ((^U- " ^^(«)) 

Using Propositions 15. ll and Proposition 15.21 it is clear that 
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along the parametrization latl = fiv'^ . Combining this with Proposition 15.31 we have that 



(5.34) ||M^(a)V^o|| < 2CJ^C^m\dX\ I - + ±.1^1^^ + 

\ 7 /if V vr fiv 



e « . 



A similar bound applies to both ||M;Vj^^(q)'(/'o|| and ||M;fc(a)-i/;o||- This completes the proof of the 
Lemma l3. II □ 
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